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Bottom Quark Mass Determination from low-n Sum Rules 

G. Corcella and A.H. Hoang ^ 

^ Max-Planck-Institut fiir Physik, Werner- Heisenberg-Institut, 
Fohringer Ring 6, D-80805 Miinchen, Germany 

We study the uncertainties in the MS bottom quark mass determination using relativistic sum rules to 0(0:3). 
We include charm mass efltects and secondary bb production and treat the experimental continuum region more 
conservatively than previous analyses. The PDG treatment of the region between the resonances T(45') and 
T(55') is reconsidered. Our final result reads: mi,{rht) 

For the sake of reliable measurements at the 
current i?-factory experiments, a precise knowl- 



edge of the bottom quark mass rUb will be essen- 
tial. In particular, the precision on the extrac- 
tion of the Cabibbo-Kobayashi-Maskawa matrix 
elements Vib from the data will depend on the 
uncertainty on the bottom mass. For example, 
an error of 60 MeV in rub leads to a 3% uncer- 
tainty in Vub from the semileptonic partial width 

In this talk we present the results of a detailed 
compilation of uncertainties in the MS bottom 
quark mass 0. Our analysis is more conserva- 
tive than an earlier one in Ref. [Sj and includes a 
number of effects that were previously neglected. 
Our method consists of determining the b mass by 
fitting the experimental moments of the bb cross 
section in e^e^ annihilation to their correspon- 
dent theoretical expressions. The moments are 
defined as follows: 



Pn = 



where R. 
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). The virtual Z contribution is strongly 
suppressed and neglected. 

We use low-n ("relativistic") moments, i.e. n < 
4. Relativistic moments exhibit the nice feature 
that they are dominated by scales of order nib 
and that they can be computed in fixed-order per- 
turbation theory. However, they have the disad- 
vantage that they strongly depend on the badly 
known experimental continuum region. As for the 
&-mass definition, we adopt the MS mass, since it 
is an appropriate definition for processes where b 
quarks are off-shell. 



(4.20 ± 0.09) GeV. 

A compilation of the theoretical moments up 
to 0(0^) was given in In our work |2] we 
also included the effects at 0{ag) of the non-zero 
charm mass and of secondary bb production, with 
the bb pair coming from gluon radiation off light 
quarks. Referring, for semplicity, to fhb(jhb) the 
theoretical moments have a simple form: 



Pn = 
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In Eq. r — rric/mb^ fi is the renormalization 
scale and we have included the contribution from 
the dimension four gluon condensate j4j . The gen- 
eral expression of P„ for fhb{fJ.) can be found in 
Refs. |2I3| . Results for the coefficients /„ are re- 
ported in Table ^ We point out that the con- 
tribution of secondary bb production and charm 
mass affects only /^°. Such effects turn out to 
be small, as can be seen comparing the values for 
/2" in Table □ with the ones of Ref. 3 . Table H 
displays the impact of charm mass corrections in 
terms of A/„ = f^°{r) - /^°(0). The smallness 
of c-mass effects is however strongly related to 
the use of the MS mass definition which we have 
adopted. In fact, if we had chosen the pole scheme 
for the bottom mass, the inclusion of the charm 
mass would have had a much stronger impact, 
as shown in Table 13 This can be understood 
from the fact that the finite charm mass repre- 
sents an infrared cut-off in the loop integrations 
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Table 1 

Coefficients of tlie theoretical moments in Eq. 



n 


1 


2 


3 


4 


A) 
J n 


0.2667 


0.1143 


0.0677 


0.0462 


riU 

J 71 


0.6387 


0.2774 


0.1298 


0.0508 


rn 

J n 


0.5333 


0.4571 


0.4063 


0.3694 




0.9446 


0.8113 


0.5172 


0.3052 


J n 


0.8606 


1.2700 


1.1450 


0.8682 


r'2'2 
In 


0.0222 


0.4762 


0.8296 


1.1240 



Table 2 

Corrections due to the non-zero charm quark mass in the MS scheme. 



r 


0.1 


0.2 


0.3 


0.4 


0.5 


A/i 


-0.0021 


-0.0078 


-0.0164 


-0.0266 


-0.0382 


A/2 


-0.0028 


-0.0091 


-0.0187 


-0.0302 


-0.0430 


A/3 


-0.0024 


-0.0101 


-0.0204 


-0.0330 


-0.0466 


A/4 


-0.0030 


-0.0109 


-0.0219 


-0.0348 


-0.0491 



and that the pole-mass definition is much more 
sensitive to infrared momenta. To evaluate the 
experimental moments, we consider the region of 
the resonances T(IS') — T(6S') and the continuum. 
We compute the moments of a generic resonance 
k in the narrow width approximation, i.e. 

{Pn)k = "TTTTTTTT^ STT+T' (^) 
a{lO GeVjmj, 

where T^.^*^ is the partial e+e~ width for the 
fc-th resonance. For the T(IS'), T(25), T(3S') 
and T{6S) we use the averages for masses and 
widths quoted in the PDG fS^. For the region 
between the T(4S') and the T(5S'), i.e. between 
10.5 and 10.95 GeV, we do not use the PDG av- 
erages, which were based on results from CUSB 
[H] and CLEO [J] Collaborations. Both exper- 
iments observed an enhancement at about 10.7 
GeV. While CUSB did not assign the enhance- 
ment to any resonance, CLEO fitted it to a " B*" 
resonance with mass ms* — 10.684 ± 0.013 GeV 
and e+e- width r^>" = 0.20 ± 0.11 keV. The 
PDG averages, on the other hand, ignore the B* 
results and, therefore, lead to a contribution to P„ 
that is smaller than the original CUSB and CLEO 
data. In our analysis we took the averages from 



the original CUSB and CLEO data, assuming the 
larger uncertainties from CLEO (See Ref. |2j for 
more details). 

As far as the continuum is concerned, we sub- 
divide it into three parts: 11.1-12.0 GeV (region 
1), where possible data may come from CLEO; 12 
GeV - AIz (region 2) and above Mz (region 3). 
There is no direct experimental data in the re- 
gion above 11.1 GeV. Nevertheless the measure- 
ments of Rb by LEP I and LEP II agree with 
the perturbative QCD prediction within 1% at 
Mz and 10% in the region between 133 and 207 
GeV explored by LEP II. It is therefore not unrea- 
sonable to rely on perturbative QCD to estimate 
the contribution to the experimental moments 
above the T(6S'). In the analysis of Ref. the 
small theoretical errors in the continuum region 
were inherently taken as the experimental uncer- 
tainties. Since this leads to an implicit model- 
dependence, we adopt a more transparent treat- 
ment and take an assigned fraction of the theo- 
retical prediction as the experimental uncertainty 
of the continuum. In this way the impact of 
the unknown experimental continuum contribu- 
tion can be traced more easily. The experimental 
moments are quoted in Tabled The uncertain- 
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Table 3 

As in Table 121 but in the pole-mass scheme. 



r 


0.1 


0.2 


0.3 


0.4 


0.5 


A/i 


0.0809 


0.1505 


0.2113 


0.2656 


0.3145 


A/2 


0.0684 


0.1267 


0.1765 


0.2203 


0.2593 


A/3 


0.0608 


0.1106 


0.1531 


0.1896 


0.2221 


A/4 


0.0545 


0.0988 


0.1358 


0.1676 


0.1764 



Table 4 

Individual contributions to the experimental moments including uncertainties. In the continuum the 
displayed uncertainties are the theoretical ones only. 



contribution 


Pi 

X 10^ GeV^ 


P2 

X 10^ GeV* 


X 10^ GeV^ 


Pi 

X 10^ GeV« 


T(IS') 


0.766(29) 


0.856(32) 


0.956(36) 


1.068(40) 


r{2S) 


0.254(16) 


0.252(16) 


0.251(15) 


0.250(15) 


T(3S') 


0.211(29) 


0.196(27) 


0.183(26) 


0.171(24) 


[T(4S) -T(5S)] 


0.251(95) 


0.218(82) 


0.190(72) 


0.165(62) 


T(6S) 


0.048(11) 


0.039(9) 


0.032(7) 


0.027(6) 


11.1 GeV- 12.0 GeV 


0.418(57) 


0.314(44) 


0.236(34) 


0.178(27) 


12.0 GeV - Mz 


2.467(26) 


0.886(21) 


0.414(13) 


0.217(8) 


Mz — oo 


0.047(1) 


0.000(0) 


0.000(0) 


0.000(0) 



ties displayed for the three continuum regions are 
the ones from the theoretical uncertainties only. 

For the determination of fhb{mb) and of the 
corresponding uncertainties, we use four meth- 
ods: we fit single moments and get directly 
mbirhb) (method 1); we determine fhb{fi) from 
single- moment fits and subsequently fhb{rhb) us- 
ing renormalization group equations (method 2); 
we fit the ratio Pn/Pn+i and get mb{fhf,) (method 
3); we determine mi,{^J.) by fitting P„/Pn+i and 
compute rhbirhb) using renormalization group 
equations (method 4). We employ four- loop 
renormalization group equations, vary the renor- 
malization scale /i between 2.5 and 10 GeV and 
use as{Mz) = 0.118±0.003, = 1.3±0.2 GeV, 
{^G^) = (0.024 ±0.024) GeV* as theoretical in- 
put. The results of our analysis are shown in Ta- 
bleini The central values for niblfhb) according to 
the four methods agree within 15 MeV. For the 
errors coming from the experimental-continuum 
regions 2 and 3 we quote both, the one coming 
from the theory uncertainties shown in Table ^ 



and the error corresponding to a 10% variation of 
the theoretical prediction. The latter error scales 
roughly linearly, i.e. assuming a 5% (20%) frac- 
tion decreases (increases) the error by a factor of 
two. In order to get the combined errors, in the 
resonance region we treat half of the errors as un- 
correlated (added linearly) and half of the errors 
as correlated (added quadratically). The errors in 
the continuum do not have any statistical corre- 
lation, hence we add them linearly. Moreover, we 
add linearly the errors coming from the resonance 
and from the continuum regions. 

We note that the errors yielded by fits of the 
first two moments Pi and P2 are rather large. 
As for the results given by fits of the moment 
ratios, the fit of P2IP?, using method 3 yields a 
rather small error of about 50 MeV. However, this 
result holds only if the same value of /i is chosen 
for both P2 and P3; a larger error would instead 
be found using independent values of /i for the 
numerator and denominator of the ratio. Since 
we believe that P3 can be calculated reliably using 
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Table 5 

Central values and uncertainties for mi,{fhh). 





Method 1 (2) 


Method 3 (4) 


n 


1 


2 


3 


1 


2 


central 


4210(4214) 


4200(4205) 


4197(4200) 


4191(4195) 


4191(4191) 


T(IS') 


14 (13) 


12 (12) 


11 (11) 


11 (11) 


9(9) 


T(2S') 


7(7) 


6 (6) 


5 (5) 


4(4) 


3(3) 


T(3S') 


14(14) 


10 (10) 


8 (8) 


7(7) 


3(3) 


45-55' 


45 (44) 


32 (32) 


22 (22) 


18 (18) 


4(4) 


T(65) 


5(5) 


3 (3) 


2 (2) 


2(2) 


0(0) 


combined 


67 (67) 


50 (50) 


38 (38) 


33 (33) 


15 (15) 


[region l]th 


27th (26th) 


17th (17th) 


Uth (Uth) 


7th (7th) 


2th (2th) 


[region 2]th 


12th(12th) 


8th (8th) 


4th (4th) 


4th (4th) 


4th (4th) 


[region 2]io% 


115 (114) 


33 (33) 


13 (13) 


49 (49) 


29 (29) 


[region 3]th 


1th (1th) 


0th (0th) 


0th (0th) 


1th (1th) 


0th (0th) 


[region 3]io% 


2(2) 


(0) 


(0) 


2(2) 


0(0) 


Srric 


0(0) 


(0) 


(0) 


0(0) 


0(0) 


Sas(Mz) 


17 (18) 


10 (11) 


6 (6) 


3(3) 


2 (2) 


Sfi 


23 (5) 


16 (14) 


11 (27) 


15 (27) 


3 (50) 


combined 


184 (166) 


77 (75) 


41 (57) 


76 (88) 


37 (85) 


total 


251 (233) 


127 (125) 


79 (95) 


110 (121) 


51 (99) 



both methods 1 and 2, we adopt the error on P3 
as our final estimate of the uncertainty in the MS 
bottom mass determination. Rounding to units 
of 10 MeV, we obtain: 

rhbifhi,) = (4.20 ± 0.09) GeV, (4) 

assuming a 10% error for the experimental con- 
tinuum regions 2 and 3. Within the error range, 
our result is in agreement with the estimate of Pj ■ 
Our error is nonetheless larger than the 50 MeV 
of Ref. 0, which is due to the different treat- 
ment of the resonance region and to the more 
conservative choice for the experimental error in 
the continuum region. 

REFERENCES 

1. M. Battaglia et al, hep-ph/0304132 

2. G. Corcella and A.H. Hoang, Phys. Lett. B 
554 (2003) 133. 

3. J.H. Kiihn and M. Steinhauser, Nucl. Phys. 
B 619 (2001) 588; Erratum-ibid. B 640 (2002) 
415. 

4. D.J. Broadhurst et al., Phys. Lett. B 329 



(1994) 103; 

P.A. Baikov, V.A. Ilyin and V.A. Smirnov, 
Phys. At. Nucl. 56 (1993) 1527 [Yad. Fiz. 
56N11 (1993) 130]. 

5. K. Hagiwara et al. [Particle Data Group Col- 
laboration], Phys. Rev. D 66 (2002) 010001. 

6. D.M. Lovelock et al., Phys. Rev. Lett. 54 
(1985) 377. 

7. D. Besson et al. [CLEO Collaboration], Phys. 
Rev. Lett. 54 (1985) 381; 

L. Garren, CLEO note CBX-84-68 (October 
1984). 



